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ABSTRACT 

From millimeter and optical observations of the Jupiter- family comet 17P/Holmes performed soon after its 
huge outburst of October 24, 2007, we derive i4n/15N=139±26 in HCN, and i4n/15n=165±40 in CN, estab- 
lishing that HCN has the same non-terrestrial isotopic composition as CN. The same conclusion is obtained for 
the long-period comet C/1995 Ol (Hale-Bopp) after a reanalysis of previously published measurements. These 
results are compatible with HCN being the prime parent of CN in cometary atmospheres. The '^N excess rela- 
tive to the Earth atmospheric value indicates that N-bearing volatiles in the solar nebula underwent important N 
isotopic fractionation at some stage of Solar System formation. HCN molecules never isotopically equilibrated 
with the main nitrogen reservoir in the solar nebula before being incorporated in Oort-cloud and Kuiper-belt 
comets. The '^C/'-'C ratios in HCN and CN are measured to be consistent with the terrestrial value. 
Subject headings: comets: general — comets: individual (17P/Holmes) — radio lines : solar system 



1. INTRODUCTION 

Comets are made of ices, organics and minerals that record 
the chemistry of the outer regions of the primitive solar neb- 
ula where they agglomerated 4.6 Gyr ago. Compositional 
analyses of comets can provide important clues on the chem- 
ical and physical processes that occurred in the early phases 
of Solar System formation, and possibly in the natal molec- 
ular cloud that predated the formation of the solar nebula 
jEhrenfreund et al.] l2005h . In particular, isotopic ratios in 
cometary volatiles are important diagnostics of how this mat- 
ter formed, since isotopic fractionation is very sensitive to 
chemical and physical conditions. However, such measure- 
ments are rare. Strong deuterium enhancements are observed 
in H2O and HCN gases of several comet comae that are char- 
acteri stic of interstellar or protosolar chemistry at low temper - 
ature ( Ehrenfreund et al.l2005tlBockelee-Morvan et al.l2005h . 
On the other hand, questions arise of how to explain the 
non-terrestria l and uniforrn value s of the ''*N/'^N isotopic 
ratio in CN dA rpigny et alj 120031: pHutsemekers et alj 120051: 
iManfroid et"ain 2005: Je hin et all |2004 . also indicative of 
interstellar-like chemi stry, while this ratio was m easured to 
be terrestrial in HCN dJewitt et alj|1997l: IZiurvse t al. 19991), 
the presumed source of CN radicals in cometary atmospheres. 
Here we present measurements of the HC^'^N/HC^^N and 
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H^^CN/H'^CN isotopic ratios acquired in comet 17P/Holmes 
by millimeter spectroscopy, together with C''*N/C^^N and 
'^CN/'-'CN measurements from visible spectroscopy. The 
brightness of this short-period (6.9 years) comet of the Jupiter 
family unexpectedly increased (from total visual magnitude 
niy = 17 to ;«,, = 2.5) on October 24, 2007, while it was at a 
distance of 1.63 AU from the Earth and at 2.44 AU from the 
Sun (lAU Circ. 8886). This huge outburst of activity, which 
is likely related to a sudden fragmentation of the nucleus fol- 
lowed by the subsequent production of a large quantity of 
grains, offered us the opportunity to search for weak spectral 
signatures of rare isotopes using complementary techniques. 

2. THE '^N/'-^N RATIO IN HCN IN 17P/HOLMES 

We carried out observations of 17P/Holmes using the 30- 
m telescope of the Institut de Radioastronomie Millimetrique 
(IRAM) located in the Sierra Nevada (Spain). Isotopic mea- 
surements were performed on October 27-28 UT. The track- 
ing of the comet was done using orbital elements K077/06 
from JPL Horizons system. The pointing of the telescope was 
checked and updated by repeated observations of a nearby 
quasar. Sky cancellation was performed by wobbling the sec- 
ondary mirror with a throw of 3' at a rate of 0.5 Hz. Four 
receivers could be operated at the same time. The 7 = 3- 
2 rotational lines of H'^C'^^N (hereafter referred to as HCN) 
at 265.9 GHz, H'^'CN (259.0 GHz) and HC'^N (258.2 GHz) 
were measured (Table [Til. The beam diameter (half-power 
beam width) of 9.5" corresponded to 11300 km at the dis- 
tance of the comet. Observations were undertaken in good 
atmospheric conditions (3-5 mm precipitable water). The 
lines were observed both at low (1 or 2 MHz) and high (62 
kHz) spectral resolutions. Spectra are shown in FigsfTKH and 
line areas are given in Table[T] Observations of HCN, H'^CN 
and HC'^N were not entirely simultaneous (Table [T]l- How- 
ever, several strong lines (CS J = 3-2, CH3CN J = 7-6 (147 
GHz) and CH3OH J = 3-2 (145 GHz) lines) were continu- 
ously observed to monitor the comet activity. Their intensities 
decreased by 12% during the period (5.5 hours long) when 
HCN, H'^CN and HC'^N data were acquired. This variation 
was taken into account when deriving the isotopic ratios. 
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TABLE 1 
Characteristics of radio lines and molecular production rates in comet 17P/Holmes. 



Molecule 



Line 



Frequency 
(GHz) 



Date UT 
(October 2007) 



Int. time Line area Opacity 

(min) (Kkms-') 



Column density 
(lO'^molcm--) 



Production rate 
(lO^f-mols-') 



HCN 

HCN 

HCN 

H"CN 

HC'^N 



y = 3-2 

/ = 3-2 F = 2-2 
J = 3-2F = 3-3 

J = 3-2 

J = 3-2 



265.886434 
265.888516 
265.884883 
259.011798 
258.156996 



27.96-27.98 
27.96-27.98 
27.96-27.98 
28.08-28.19 
27.96-28.19 



15 
15 
15 
90 
145 



18.23 ±0.14 
0.958 ± 0.083 
1.097 ±0.073 
0.220 ± 0.037 
0.191 ±0.021 



0.56 

0.03 

0.03 

0.008 

0.006 



20.8 ±0.17 
20.7 ± 1.8 
23.7 ± 1.6 
1.74 ±0.30 
1.51 ±0.16 



19.85 ±0.16 
19.74 ± 1.71 
22.60 ±1.50 
0.164 ±0.028 
0.140 ±0.015 



Note. — The integration time is on+off source; error bars are l a; line are as refer to the main beam brightness temperature scale; column densities and 
production rates were determined using the radiative transfer model of lBiver et a l. 1 1999) with x,ie = 0.5 and a steady state coma with Ti-,,, = 45 K as described in 



components of H CN J = 3-2 are shown. 



the text. The production rate determined from the whole HCN line (resp. HCN hyperfine components, H CN and HC 
using Tj,,, = 65 K derived by Dello Russo et al. (2008) from HCN infrared observations on October 27.6, 2007. 

In contrast to the HC'^N and H'^CN lines, the HCN J = 3-2 
is optically thick (Table[T]i. Therefore, radiative transfer mod- 
eling is required to reti-ieve the HCN/H'^CN, HCN/HC'^N 
isotopic ratios if the HCN J = 3-2 line is used in the analysis. 
The coma of comet Holmes was in a non-equilibrium regime 
following its outburst, with gas phase species released by icy 
grains. Unavoidable simplifying model assumptions on the 
coma structure could introduce systematic opacity-dependent 
errors in the retrievals and affect the determination of the iso- 
topic ratios. However, the HCN J = 3-2 line is split into six 
hyperfine components, two of which (F = 2-2, F = 3-3) are 
well separated from the core of the line, and detected in the 
HCN spectrum (Fig.|2] Table[T]). These hyperfine components 
have intrinsic line strengths of 3.7% the total strength accord- 
ing to hyperfine statistical weights, and are optically thin (Ta- 
ble [B- HCN/H'^CN, HCN/HC'^N abundance ratios deter- 
mined using the HCN optically thin hyperfine components do 
not depend upon model assumptions on coma temperature, 
structure and temporal variability, because in this case emis- 
sion lines from molecules in the same excitation state and in 
the same regions of the coma are compared. They are di- 
rectly given by the line intensity ratios, albeit minor correc- 
tions accounting for slightly different line frequencies. Using 
the HCN hyperfine components, and correcting for the non- 
simultaneity of the HCN, H'^CN and HC'^N measurements, 
we derive H'^CN/H'^CN = 114 ± 26 and HC^^N/jjciSp^ ^ 
139 ±26. 

Column densities and pr oduction rates dete rmined with the 
radiative transfer model of lBiver et alJ (1 19991) are given in Ta- 
ble [U We assumed a steady-state isotropic parent molecule 
distribution, with a gas velocity of 0.56 km s"', and a gas 
kinetic temperature Tki„ of 45 K inferred from IRAM obser- 
vations of mutiple lines of CH3OH on October 29.0 UT. In- 
terestingly, similar isotopic ratios are obtained using the HCN 
production rate deduced from the whole J = 3-2 line as with 
the model independent hyperfine lines method. As shown in 
Fig.|2l when the day and night side HCN velocities are fixed 
to 0.6 and 0.4 km s"', the model provides a satisfactory fit 
to the shape of the HCN line. This suggests that our descrip- 
tion of the HCN spatial distribution and excitation is correct in 
first approximation. This conclusion is supported by the good 
agreement between the HCN production rate measured on 
Octob er 27.6 in a smaller (~ 1") aperture (iDello Russo et alJ 
I2OO8I) and those reported in this work. 

3. THE '*N/'5n RATIO IN CN IN 17P/HOLMES 

High-resolution optical observations were performed to 
measure the '**N/'^N and '^C/'-'C ratios in CN. Spectra of 
the B2I]+-X2I]+(0,0) CN band at 388 nm were obtained on 
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Fig. 1 .— Spectra of the J = 3-2 lines of HC'^^N and H'^CN in comet 
17P/Holmes on 27-28 October 2007. The velocity frame is with respect to 
the comet rest velocity. The positions and relative intensities of the hyperfine 

tI3 




Fig. 2. — Model fitting (continuous line, see text) to the J = 3-2 HCN 
line profile observed in comet 17P/Holmes on Oct. 27.97 UT (dotted line). 
Positions and relative intensities of the hyperfine components are shown. 

October 25.4, 28.3, 29.4, 30.4, 31.4 and November 18.4 and 
19.3 2007 UT with the 2DCoude spectrograph at the 2.7-m 
Harlan J. Smith telescope of the McDonald Observatory. A 
series of short exposures, from 30 seconds to 5 minutes, were 
also collected on October 29.6 2007 UT with the High Reso- 
lution Echelle Spectrometer (HIRES) of the Keck 1 telescope 
installed on Mauna Kea (Hawaii). The observations were car- 
ried out under clear weather and low airmass. The slit of 
the spectrograph was in both cases ^1" wide and ~7" long 
providing a resolving power of about A/AA = 60 000 (0.03 
A/pixel). The slit was centered on the false nucleus in the 
case of the Keck spectra and displaced in the coma (by up to 
20") for the McDonald exposures (20 min) to reduce the con- 
tamination by the strong dust-reflected spectrum. The dust- 
reflected sunlight underlying the spectral lines (among which 
are ro-vibrational Hnes of '"C'^^N, '^^C'-^N and '^c'^N (0,0) 
band) was removed by subtracting a solar reference spectrum 
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Fig. 3. — A small part of the co-added Keck and McDonald CN (0,0) band 
spectra of comet 17P/Holmes showing seven R branch lines (R3 to R9) of CN 
isotopes (black line). The dust-scattered solar spectrum was removed. The 
positions of the "C'''N (red ticks) and '-C'^N (blue ticks) lines are indicated; 
the con'esponding lines of the main isotope are on the left. The synthetic 
spectrum with '^C/"C=90 and '■*N/'''N = 165 is superimposed (blue line). 

after the appropriate Doppler shift, profile fitting and normal - 
ization were applied (lArpignv et al.ll2003t iJehin et al]|2004t) . 
The individual CN (0,0) spectra were then combined with an 
optimal weighting scheme in order to maximize the overall 
signal-to-noise ratio. Synthetic spectra of '^C''*N, '■'C''*N, 
'^C'^N were compute d for each observing circu mstance using 
a fluorescence model (IZucconi & Festoulll985l) . We took into 
account slightly different excitation conditions over the period 
of the observations caused by small variations of the heliocen- 
tric distance and velocity. Collisional e ffects were empiricall y 
estimated by fitting the '^C''*N fines (iManfroid et al.ll2005h . 
The synthetic spectra were then co-added in the same way as 
the data. The isotope mixture was adjusted to best fit the ob- 
served continuum-subtracted spectrum. We considered seven 
R branch lines (R3 to R9), as shown in Fig. [3] The isotopic 
ratios of '^C/'^C and '''N/'-''N in the final co-added Keck and 
McDonald spectrum are estimated to be 90 ± 20 and 165 ± 
40, respectively. 

4. REANALYSIS OF THE COMET HALE-BOPP DATA 

The '"^N/'^N ratios measured in HCN and CN in comet 
17P/Holmes are consistent with each other In contrast, an 
HC''*N/HC'^N ratio marginally higher than the Earth atmo- 
spheric value (272) w as reported for comet C/1995 01 (Hale- 
Bopp) (iJewittetal .1 11997: Ziurysetal. 1999), to be com- 

S iared to the value '"^N/'^N = 140± 35 measured in CN 
Arpignv et al.ll2003h . This led us to reanalyse the measure- 
ments made in comet Hale-Bopp; 

Our reanalysis of the data of IZiurys et al.l (1 19991) obtained 
on March 24 and 25, 1997 yields '^C/'^C = 65 ± 13 and 
i4n/15]^ = 152 ± 30 (with a 10% calibration uncertainty in- 
cluded) as compared with the HC'"^N/HC'^N production rates 
ratio of 100 ± 20 and 286 ± 82 given by these authors. Ziurys 
et al. used approximate formulas to analyse the optically thick 
HCN lines, assuming in addition an inappropriate value for 
the rotational temperature of HCN. Instead, we used full ra- 
diative transfer modeling. Our determinations are likely more 
reliable as the HCN production rates that we derive from the 
J = 1-0, 2-1 and 3-2 HCN fines observed on March 24 are 
consistent within 5%, while the values inferred by Ziurys et 
al. from the different lines differ by up to a factor of 2. The 
HC'^N and H'^CN J = 3-2 lines were observed on March 25. 
As one cannot exclude day-to-day variations in HCN produc- 
tion from March 24 to 25, the HC'^N/H'^^CN ratio of 0.43 ± 
0. 10, deduced from March 25 data only, could be more secure. 
In this case, assuming '^C/'^'C to be equal to the terrestrial 
value of 89, we deduce '''N/'^N = 207 ±4 8. 
The data obtained bv lJewitt et"ai] d 1997b at flie James Clerk 



Maxwell Telescope (JCMT) are public and available from the 
Canadian Astronomy Data Centre. Reanalysing these data, 
we found: 1) the HCN data were acquired during night time 
near sunrise, while the H'^CN line and especially the HC'^N 
line were observed later during daytime; the beam efficiency 
may have degraded by 20% according to JCMT specifica- 
tions; 2) the HCN line was observed with the receiver B3 
tuned in double sideband (DSB), while the other lines were 
observed in single sideband (SSB); HCN spectra of calibra- 
tion sources obtained with the same receiver DSB tuning as 
used for the cometary observations show signals in excess 
of 15% with respect to reference spectra of the sources; 3) 
the H'^CN J = 4-3 line is ble nded with the S O? 13? 12- 
12i.ii line at 345.338538 GHz dLis et all [19971) ; using the 
SO2/HCN productio n rate rati o determined in comet Hale- 
Bopp (Bockelee-Mory an et al,' 2000), we estimate that it af- 
fects the H'^CN fine intensity by 20%; 4) more critical, the 
HC'^N spectra present scan-to-scan intensity variations by a 
factor of 2 (a factor of 10 above the fluctuations related to the 
statistical noise) that are likely of instrumental origin. Tak- 
ing these corrections into account, we infer '^C/'- ^C = 94 ifc 8 
and^N/'^N = 205 ± 70, while the values given in lJewitt et all 
(119971) are 100 ± 12 and 323 ± 46, respectively. The large un- 
certainty in our '"^N/'^N determination reflects the dispersion 
of the HC'^'N measurements. 

We conclude that the '''N/'^N ratio in HCN is rather un- 
certain for comet Hale-Bopp but is consistent with the value 
measured in CN. 

5. IMPLICATIONS 

The '^C/'^'C values in HCN and CN are in agreement 
with the terres trial value of 89 and previou s measurements 
in comets (e.g.. lBockelee-Morvan et al]|2005l) . 

The '^N/'^N ratios measured in HCN and CN in comet 
17P/Holmes both correspond to a factor of two '^N enrich- 
ment relative to the Earth atmospheric value, and are con- 
sistent with the '''N/'^N ratios measure d in CN in a dozen 
comets which cluster a t 141 ± 29 dAjpignv et al.l l2003t 
Hutse mekers et al.1 l2005t IManfroid et al.1 l2005t iJehin et alj 
2004). The discrepancy between HCN and CN isotopic ratios 
previously found for comet Hale-Bopp led to the suggestion 
of a CN production mechanism other than HCN photolysis 
in cometary atmospheres, possibly from the thermal degra- 
dation of '^N-rich refractory organics present in dust grains. 
This interpretation was supported by the presence of CN jets, 
by the radial distribution of CN, found to be generally less ex- 
tended than expected from HCN photodissociation, and by the 
CN/HC N production rat e ratio which exceeds unity in some 
comets jFrav et al.l2005h . However, questions arose of how to 
explain the equally low C''*N/C'^N value observed in comet 
Hale-Bopp at large heliocentric distance where the CN radi- 
cals w ere expected to be rnainly HCN ph otodissociation prod- 
ucts (.Manfroid etani2005HRauer et al.ll2003.) . Our reanalysis 
of the Hale-Bopp data, which shows that the '^N/'^N ratio in 
HCN encompasses the CN value, solves this issue. Also, a 
similar isotopic ratio in HCN and CN provides a better expla- 
nation to the uniform values of the C'^N/C'^N ratio among 
comets which exhibit large differences in dust-to-gas produc- 
tion rate ratios. 

Our isotopic measurements are compatible with HCN being 
the prime parent of CN in cometary atmospheres. For comet 
Hale-Bopp, the production r ates of the two species were 
found approximately equal feRauer et al.l 120031: iFrav et all 
l2005h . The complex and variable structure of 17P/Holmes's 
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coma after its outburst makes comparisons of the HCN and 
CN production r a tes di fficult for this comet. However, 
iDello Russo et al.l (|2008') note that the good agreement be- 
tween the HCN/H2O abundance determined from infrared 
spectra and the OH/CN abundance ratio measured with nar- 
rowband photometry (Schleicher 2007) is consistent with CN 
being mainly produced by HCN photolysis. Yet, we cannot 
exclude that, in other comets, CN has other major progenitors 
(dust or gas-phase species) sharing the same low '"^N/'^N iso- 
topic ratio, which would mean that '^N enrichment is possibly 
a general property of CN-bearing compounds in comets. 

Comets Hale-Bopp and 17P/Holmes are of different dy- 
namical families, the former originating from the Oort cloud, 
and the latter from the trans-Neptunian scattered disc. CN 
anomalous nitrogen isotopic composition is observed in 
a number of comets of these two dynamical populations. 
Hence, comets issued from these two reservoirs likely exhibit 
similar anomalous N isotopic composition in HCN. 

The '^N excess measured in cometary HCN (S^^N ^ 1000 
per mil relative to the Earth atmospheric value) and CN 
daughter-product cannot result from isotopic fractionation in 
the comet atmosphere. It is comparable to the extr eme enrich- 
ments measured in interplanetary dust particles (Flo ss et al 
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2006) (IDPs) and carbonaceous meteorites (Buseman netal. 



2006). High (5'^N values are also present at the submicrome- 
tre scale in the dust par ticles collected by the S tardust mission 
in comet 8 IP/Wild 2 (iMcKeegan et alj|2006l) . In IDPs and 
meteorites, the '^N-rich nitrogen is carried by non-volatile 
macromolecular organic material and is generally believed to 
be a remnant of interstellar chemistry, though a nucleosyn- 
thesis source is often considered as an alternative mechanism 
given the presence of presolar gr ains with high '^N excesses 
in the matrice of meteorites (e.g..l ZinneiJll998h . These com- 
plex organics possibly formed from UV or cosmic radiation 
processing of simple ices in the presolar cloud or, at later 
stages, in the cold regions of the solar nebula. The enrich- 
ments reported here provide the first evidence for the pres- 
ence of high '^N anomalies in the volatiles that composed the 
icy phase of the outer solar nebula, possibly representing the 
precursors of the primitive refractory organics. 

Cometary HCN and H2O ices exhibit strong en- 
richments in deuterium with respect to the cosmic 
D/H value jEhrenfreund et al.ll2005t iBockelee-Morvan et al] 
l2005h . which are believed to reflect ion-molecule and gas- 



grain fractionation reactions that took place at low tempera- 
ture in the early phases of Solar System formation or in the 
natal molecular cloud. The interpretation of the '^N enrich- 
ment in HCN (by a factor of 3 with respect to the pr otoso- 
lar value in th e mai n nitrogen reservoir; Fouchet et al. I2004t 
Meibom et al. |2007l) is not as compelling as for deuterium, be- 
cause there is still little evi dence for N isotop ic fractionation 
in the interstellar medium (llkeda et alj|2002l) and predicted 
'^N enhancements in HCN for exchange re actions involv- 
ing th e main nitrogen reservo i r N2 a r e mod est dTerzieva et al.l 
|2000b. iRo dgers & Charnle vI (|2004 |2008|) show fliat highly 
fractionated NH3 ice could have formed in interstellar or pro- 
tosolar material if N2 was converted into atomic nitrogen. 
While this mechanism is attractive to account for the high 
'^N excesses in primitive refractory organics, if synthesized 
from NH3, the chemical network for explaining the HCN iso- 
topic anomaly in comets still has to be proposed. Alterna- 
tive mechanisms include a nucleosynthetic origin and pho- 
tochemical self-shielding in the solar nebula, similar to that 
proposed for explaining the oxygen isotope anomalies in me- 
teorites (Clavton 2002a b). Irrespective of the mechanism, 
protosolar HCN never isotopically equilibrated with nebular 
gas at later phases of Solar System evolution. 
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